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Abstract 

We investigate what would happen to the time dependence of a pulse re- 
flected by a disordered single-mode waveguide, if it is closed at one end not 
by an ordinary mirror but by a phase-conjugating mirror. We find that the 
waveguide acts like a virtual cavity with resonance frequency equal to the 
working frequency loq of the phase-conjugating mirror. The decay in time of 
the average power spectrum of the reflected pulse is delayed for frequencies 
near wq. In the presence of localization the resonance width is exp(— L//), 
with L the length of the waveguide, I the mean free path, and Tg the scattering 
time. Inside this frequency range the decay of the average power spectrum is 
delayed up to times t ~ exp(L//). 

PACS numbers: 42.65.Hw, 42.25.Dd, 72.15.Rn 
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I. INTRODUCTION 



The reflection of a wave pulse by a random medium provides insight into the dynamics 
of localization [jlj-Q. The reflected amplitude contains rapid fluctuations over a broad range 
of frequencies, with a slowly decaying envelope. The power spectrum a{uj,t) characterizes 
the decay in time t of the envelope at frequency u. In an infinitely long waveguide (with 
propagating modes), the signature of localization 

(a(cj,t)) oc for t > A^V„ (1) 

is a quadratic decay of the disorder-averaged power spectrum (a), that sets in after A^^ 
scattering times r^. 

The decay (0) still holds over a broad range of times if the length L of the waveguide is 
finite, but much greater than the localization length ^ = (A^ + 1)/ (with / = cTg the mean 
free path). What changes is that for exponentially large times t ^ Tg exp{L/l) the quadratic 
decay becomes more rapid oc exp(— constant x In^t). This is the celebrated log-normal tail 
IfzHTli. We may assume that the finite length of the waveguide is realized by terminating 
one end by a perfectly reflecting mirror, so that the total reflected power is unchanged. 

In this paper we ask the question what happens if instead of such a normal mirror 
one would use a phase-conjugating mirror [^,0. The interplay of multiple scattering by 
disorder and optical phase conjugation is a rich problem even in the static case ||I^-[16|. 
Here we show that the dynamical aspects are particularly striking. Basically, the disordered 
waveguide is turned into a virtual cavity with a resonance frequency uq set by the phase- 
conjugating mirror. 

We present a detailed analytical and numerical calculation for the single-mode case (A^ = 
1). For times t 3> we find that a{uj,t) has decayed almost completely except in a narrow 
frequency range oc exp{—L/l) around uq. In this frequency range the decay is delayed 
up to times t c::^ Tg exp(L//), after which a log-normal decay sets in. The exponentially large 
difference in time scales for the decay near ujq and away from ujq is a signature of localization. 



II. FORMULATION OF THE PROBLEM 

A. Scattering theory 

A scattering matrix formulation of the problem of combined elastic scattering by disor- 
der and inelastic scattering by a phase-conjugating mirror was developed by Paasschens et 



al. |T^. We summarize the basic equations for the case of a single propagating mode in the 
geometry shown in Fig. ^ A single-mode waveguide is closed at one end {x = 0) by either a 
normal mirror or by a phase-conjugating mirror. Elastic scattering in the waveguide is due 
to random disorder in the region < a; < L. For simplicity we consider a single polarization, 
so that we can use a scalar wave equation. 

The phase-conjugating mirror is pumped at frequency ujq. This means that a wave 
incident at frequency ujq + u will be retro- reflected at frequency uo — uj, for u <^ uq. For 
X ^ L the wave amplitude at frequencies uj± = uq ± u is an incoming or outgoing plane 
wave, 
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FIG. 1. The geometry under investigation consists of a single-mode waveguide with a mirror 
at X = 0. It can be a normal mirror or a phase-conjugating mirror. There are randomly positioned 
obstacles between x = and x = L. 

M±(r, t) = Re (j)^^ exp[—ik±{x — L) — iio±t\il)±{y, z), (2a) 
u°^\r, t) = Re (f)"f exp[ik±{x - L) - io^±#±(|/, z). (2b) 

Here k± = fco ±cj/c is the wavenumber at frequency uj±, with the wavenumber at uoq and 
c = du/dk the group velocity. The transverse wave profile '4'±{y, z) is normalized such that 
the wave carries unit flux. 

The reflection matrix relates the incoming and outgoing wave amplitudes, according to 



out / \ / I \ iri 

(3) 



\<P- J V ''^ 

The reflection coefficients are complex numbers that depend on u. They satisfy the symme- 
try relations 

rl_(u;) = r++(-a;), r%{u) = r+_(-a;). (4) 
If there is only reflection at the mirror, and no disorder, then one has simply 



r / I — e" 



2ik-L (5) 



for a normal mirror, and 



for a phase-conjugating mirror operating in the regime of ideal retro- reflection. (We will 
assume this regime in what follows.) 

We wish to determine how the reflection coefficients are modified by the elastic scattering 
by the disorder. For this we need the elastic scattering matrix 

(;::)• « 

The reflection coefficients r,r' and transmission coefficients t,t' describe reflection and trans- 
mission from the left or from the right of a segment of a disordered waveguide of length 
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L. The matrix 5* is unitary and symmetric (hence t = t'). The basis for S is chosen such 
that r = r' = 0, t{±uj) = e^^^^ in the absence of disorder. The relationship between the 
coefficients in Eqs. (|^) and (|^ is 

r++(cj) = r\uj) + t{uj)[l - r*{-uj)r{uj)]-\*{-uj)t{uj), (8a) 
r+_{uj) = -it{uj)[l - r*{-uj)riuj)]-h*{-uj), (8b) 

for a phase-conjugating mirror. For a normal mirror there is no mixing of frequencies and 
one has simply 

r++(cj) = r'(cj) - t{uj)[l + r{u;)]-H{uj), (9a) 
r+_(cj) = 0. (9b) 

In each case the matrix of reflection coefficients is unitary, so 

|r++(a;)|2 + |r+_(cu)|2 = l. (10) 



B. Power spectrum 



We assume that a pulse oc 6{t) is incident at x = L [corresponding to 0™ = 1 for all lj 



in Eq. (|^)]. The reflected wave aX x = L has amplitude 

Mout(t) = Re e-^"«* 



°^ da; 
2^ 



r__(uo)+r*{uo] 



(We have suppressed the transverse coordinates y,z for simplicity of notation.) Using the 
symmetry relations (^, we can rewrite this as 



Uoutit) = Re e' 



-iujot 



l-oo 27r 

The time correlator of the reflected wave becomes 



-iujt 



(12) 



«out(t)nout(t + t') = \ Re e*-«*' I" ^ |_ 



X 



, do. do/ ^^^^.^^-^t^iu't' 
oo 27r 



(13) 



plus terms that oscillate on a timescale 1/ujq. We make the rotating wave approximation 
and neglect these rapidly oscillating terms. The power spectrum a of the reflected wave is 
obtained by a Fourier transform, 



a{(jj, t) 



oo 



dt' cos [(u;o + Ld)t']Uont{t)Uout{t + t') 



Re 



iSujt 



-oo 27r 

where we have introduced the correlator in the frequency domain 



(14) 



4 



a{u!, Su) 



Integration of the power spectrum over time yields, using also Eq. (pUj), 

/oo 
dt a{uj, t) = Re a(u;, (5ti; = 0) = | + | Re r+ 
-oo 



(15) 



(16) 



For a normal mirror r+_(co') = and a{uj, ^cu = 0) = |, expressing flux conservation. For the 
phase-conjugating mirror there is inelastic scattering, which mixes the frequency components 
uj and —uj. The constraint of flux conservation then takes the form 



a{uj, 5uj = {)) + a(— cj, 5uj = 0) 



(17) 



This follows from the symmetry relations (|^) and the unitarity of the reflection matrix. 
Eq. ( p!7D implies that a{uj = 0,6uj = 0) = |. 



III. RANDOM SCATTERERS 

We assume weak disorder, meaning that the mean free path I is much larger than the 
wavelength 27r/fco- The multiple scattering by disorder localizes the wave with localization 
length equal to 21. We consider separately the case of a phase-conjugating mirror and of a 
normal mirror. 



A. Phase-conjugating mirror 



We flrst concentrate on the degenerate regime of small frequency shift u, and will simplify 
the expressions by putting u = from the start. We note that r++(0) = 0, 'r+_(0) = —i, as 
follows from Eq. (|^) and unitarity of the scattering matrix (|^). Using Eqs. and (0), we 
arrive at the power spectrum in the frequency domain 



a{0,6uj) = ^-(r'{6uj) + [l - r*{-6uj)r{6uj)] ^ [t\6ujy {-6uj) - it{Suj)t* {-6uj) 



;i8) 



The scattering amplitudes have the polar decomposition r = V Rexp{i9), r' = V Rexp{i6'), 
t = — Rexp[^i{9 + 9')], with R, 9, 9' real functions of frequency. The phase 9' may be 
assumed to be statistically independent of i?(±5cij), 6'(±(5co'), and uniformly distributed in 
(0,27r). (This is the Wigner conjecture, proven for chaotic scattering in Ref. Ijl^-) In this 
way only the last term in Eq. (|T^ survives the disorder average (■■■), 



4(a(0,(5c<;) 



t{5uj)t*{-5uj) 
1 — r*{—5uj)r{5uj) 



HZ. 



(19) 



m=0 



where we have deflned Zr, 



{t{5u)t*{-5uj)[r*{-5u)r{5uj)Y'). 



The moments Zm satisfy the Berezinskii recursion relation [p!8l , |l9 



^dZm 



m'^{Zm+i + Z, 



m— 1 



2Zm) + (2m + l){Z„i+i — Zm) + liTsbioi^m + 



(20) 
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with Ts = l/c the scattering time. [The mean free path I accounts only for backscattering, 
so that the scattering time in a kinetic equation would equal ^t^.] The initial condition 
is Zm{L = 0) = Sm,o- In App. A we derive an analytical result for {a{0,6uj)) in the small 
frequency range \n{l/TsSuj) > L/l ^ 1. It reads 

(a(0,M) = i dA; ik {-2iTs5ujy''-^/^2-^"'-^/^T^{l + ik)T{l - ik) 

X r~\l + ik)r~\tk) exp[-(i + k^)L/l]. (21) 

The initial decay is determined by the contributions of the poles at k = —\h ~\h 

a(0,(5cj)) = \ + \iTjujexp{2L/l) - ^t^6uj^ exp{QL/l) + 0{6uj^). (22) 

The result (|2lD is plotted in Fig. ^ for L/l = 12.3. We compare with the data from a 
numerical solution of the wave equation on a two-dimensional lattice, using the method of 



recursive Green functions [20|. (The method of simulation is the same as in Ref. [15|, and 
we refer to that paper for a more detailed description.) The agreement with the analytical 
curves is quite good, without any adjustable parameter. The 5u;-dependence of {a{0,6uj)) 
for large L/l occurs on an exponentially small scale, within the range of validity of Eq. (pH). 
A Fourier transform of Eq. (p]) yields the average power spectrum in the time domain 



for \n{t/Ts) > L/l > 1, with the result 

[a{0, t)) = l^'/^L/l)-^/^ exp{-L/4l)T;'/H-^/^ ln(4t/r,) exp [ - (l/AL) \n''{4t/T,)] . (23) 



The leading logarithmic asymptote of the decay is log- normal oc exp[— (//4L) In^t], charac- 
teristic of anomalously localized states |]7HTT|. 

These results are calculated for = 0, and remain valid as long as ^ r^^ exp(— L//). 
This is the degenerate regime. For larger frequency mismatch u one enters the non- 
degenerate regime. The power spectrum in that regime is the same as for a normal mirror, 
calculated in the next subsection. 



B. Normal mirror 

For comparison we discuss the known results for a disordered waveguide connected to a 
normal mirror instead of a phase- conjugating mirror. Since r+_ = 0, one has from Eq. (|15D 

4(a(cj, Suj)) = (r++{uj + Suj)rl^{iu)) = Ri. (24) 
The quantities Rm = ([r++{uj + 6uj)r^^{uj)]"^) satisfy the Berezinskii recursion relation 



dR 

m+l 

+ Rm-i — 2Rm) + 2iTs6ujmRm- (25) 
The initial condition is Rm{L = 0) = 1 for all m. The solution for \n{l/Ts5uj) ^ L/l is 



known Ell], and gives the average power spectrum 
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FIG. 2. Average power spectrum for reflection by a disordered waveguide {L/l = 12.3) con- 
nected to a phase-conjugating mirror [solid curves, from Eq. (|2l])] or a normal mirror [dashed curves, 
from Eq. (28)]. The data points follow from a numerical simulation. There is no adjustable param- 
eter in the comparison. Notice the much faster frequency dependence for the phase-conjugating 
mirror (top panels), compared to the normal mirror (bottom panels). 



^a{uj, Suj) 
1 

+ - 



-2iTsSuj[ Ki 



IT J-oo 



dA; ksinh{7rk){l + k^y^K 



2ik 



exp 



i + k')L/l 



(26) 



with K a Bessel function. [The result ( PBD does not require L/l ^ 1, in contrast to Eq. (PT|).] 
The initial decay is dominated by the contributions of the poles at = — —^h 



a{uj, 6uj) 



\iTs5ujL/l - ir25cj2exp(2L//) + 0(5cu^ 



(27) 



Comparison of Eqs. ( pq ) and ( PTf ) with Eqs. (^) and ( p^ ) shows that the decay is much 
slower for a normal mirror than for a phase-conjugating mirror. The characteristic frequency 
scale is larger by a factor exp(2L/Z). So Eq. (|26D is not sufficient to describe the entire decay 
of {a{(jj,5uj)), which occurs in the range Ts5uj ^ 1. The decay in this range is obtained by 
putting the left-hand-side of Eq. (p5[) equal to zero, leading to |5|,p^ 



a(u;, 5uj) 



(28) 
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where Ei is the exponential integral function. The range of validity of Eq. (pS]) is 
\n{l/Ts6uj) <C L/l and L/l ^ 1. The result (|28|) is plotted in Fig. ||, and is seen to agree 
well with data from the numerical simulation. 

For ln(t/rs) <^ L/l (and L/l S> 1) one can perform the Fourier transform of Eq. (^) to 
get the average power spectrum in the time domain 

{a{uj, t)) = lTs{t + 2r,)-2, t > 0. (29) 

It decays quadratically oc for t/vg ^ 1. For exponentially long times, t ^ Ts exp(L//), 
one should instead perform the Fourier transform of Eq. (|26|). One finds that the quadratic 
decay crosses over to a log- normal decay oc exp[— (//4L) In^ t] 0], the same as for the phase- 
conjugating mirror. 



IV. CONCLUSION 

We have shown that the interplay of phase-conjugation and multiple scattering by dis- 
order leads to a drastic slowing down of the decay in time t of the average power spectrum 
{a{uj,t)) of frequency components u of a reflected pulse. The slowing down exists in a nar- 
row frequency range around the characteristic frequency uoo of the phase- conjugating mirror 
(degenerate regime). If u is outside this frequency range (non-degenerate regime), the power 
spectrum decays as rapidly as for a normal mirror. 

The slowing down can be interpreted in terms of a long-lived resonance at ujq, that is 
induced in the random medium by the phase-conjugating mirror. This resonance is known 
from investigations of the static scattering properties |]TB|. The resonance is exponentially 
narrow oc r~^exp(— L/Z) in the presence of localization (with the scattering time, L the 
length of the disordered region, and / the mean free path). The resonance leads to the 
exponentially large differences in time scales for the decay of the power spectrum in the 
degenerate regime and the non-degenerate regime. 

We have restricted the calculation in this paper to the case of a single propagating mode, 
when a complete analytical theory could be provided. We expect that the A^-mode case is 
qualitatively similar: An exponentially large difference in time scales oc exp(L/^) for the 
decay in the degenerate and non-degenerate regimes provided the medium is localized [L 
large compared to the localization length C, = {N + 1)/]. In the diffusive regime we expect 
(a(a;,t)) to decay on the time scale of the diffusion time Ts{L/1)^. The difference with the 
non-degenerate regime (or a normal mirror) is then a factor {L/iy instead of exponentially 
large. 

In final analysis we see that phase conjugation greatly magnifies the difference in the 
dynamics with and without localization. Indeed, if there is no phase-conjugating mirror 
the main difference is a decay oc in the diffusive regime versus in the localized 

regime 0, but the characteristic time scale remains the same (set by the scattering time 
Ts). We therefore suggest that phase conjugation might be a promising tool in the ongoing 



experimental search for dynamical features of localization |23,24 
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APPENDIX A: POWER SPECTRUM IN THE FREQUENCY DOMAIN 



We show how to arrive at the resuh (PTI) starting from the recursion relation (PO). We 
assume hi{l/Ts6uj) > L// ^ 1. It is convenient to work with the Laplace transform 



ZM)= exp(-AL/Z)Z„(L) 

JO I 



(Al) 



of the moments Zm- The recursion relation (1201) transforms into 



XZm{\) - Sm,0 = Zm+l{\) + ^m-l(A) - 2Zm{\) + (2m + 1) Z.m+l{\) - ^m(A) 

-/3(2m+l)Z^(A), 

with P = —2iTs5uj. 

For small \(3\ and large m this equation can be written as a differential equation, 



(A2) 



2(92Z(m,A) , ^_dZ{m,\) 



+ 2m- 



{\ + 2pm)Z{m,X) = 0, 



dm? dm 

where m is now considered to be a continuous variable. The solution of Eq. { \A3\ ) is 

Z{m,X) = C{X,P){Pm)-'/'K^^ (^2^2^^ 



(A3) 



(A4) 



The factor C(A, /?) is determined by matching to the solution of Eq. ([A2| ) for l3m 0, 
m — > cxo, that has been calculated in Ref. [^. The result is 



C(A,/3) = Ait(3^/^t[1 + lVmx)T-'{l + lVmx)T-' [^VT+AX 
X exp fiVl + 4Aln(/3/8) 



(A5) 



To obtain the power spectrum ([I9D we replace the sum over m by an integration, with 
the result 



oo 

Y: zux) = 2'/\p-'/'T\l + iv/rT4A)r(i - iv^rT4A)r-i(i + IVi + AX 



m=0 



X r-MiVl + 4A)exp |Vl + 4Aln(^/8) 



(A6) 



There are poles at A = n{n + 1), n = 0, 1, 2, . . . and a branch cut starting at A = —1/4. 
When doing the inverse Laplace transform we put the integration path in between the poles 
and the branch cut. The final result is given by Eq. (^). The reason that we need the 
condition L/l ^ 1 is that Eqs. ( [A4| ) and ( [^ ) are only correct for m ^ 1. The first terms 
in the sum J2m=o are important for L/l ^ 1, but can be neglected for L// ^ 1. 
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